Observation of Weak-Limit Quasiparticle Scattering via Broadband Microwave 
Spectroscopy of a d-Wave Superconductor 



(N 
O 

o 
< 



S3 
O 
O 

Or 
c3 



I 

C 

O 

o 



> 

m 
in 
m 



I 

o 
o 



P. J. Turner, R. Harris, Saeid Kamal, M. E. Hayden,| D. M. Broun, D. C. Morgan, A. 
Hosseini, P. Dosanjh, G. Mullins, J. S. Preston,[j] Ruixing Liang, D. A. Bonn, and W. N. Hardy 
Department of Physics and Astronomy, University of British Columbia, Vancouver, B.C., Canada V6T 1Z1 

(Dated: February 1, 2008) 

There has long been a discrepancy between microwave conductivity measurements in high temper- 
ature superconductors and the conductivity spectrum expected in the simplest models for impurity 
scattering in a d-wave superconductor. Here we present a new type of broadband measurement 
of microwave surface resistance that finally shows some of the spectral features expected for a 
d x 2_ y 2 pairing state. Cusp-shaped conductivity spectra, consistent with weak impurity scattering 
of nodal quasiparticles, were obtained in the 0.6-21 GHz frequency range in highly ordered crystals 
of YBa2Cu3O6.50 and YBa2Cu3 06.99- 

PACS numbers: 74.25.Nf, 74.72.Bk 



It is now widely accepted that cuprates exhibiting high 
temperature superconductivity involve a pairing state 
with d x 2_ y 2 symmetry Q. However, there is considerable 
debate over the extent to which conventional Bardeen- 
Cooper-Schrieffer theory can account for the properties 
of this superconducting state, a debate that is closely 
related to the difficulty in understanding the peculiar 
metallic state in these materials. The question is, de- 
spite the new physics being encountered in the cuprate 
phase diagram, is the d x 2_ y 2 superconducting state sim- 
ple? In this letter we present the first low frequency 
conductivity spectra that answer this in the affirmative, 
showing an evolution at low temperatures towards the 
cusp-shaped spectrum expected for weak-limit impurity 
scattering in a clean d-wave superconductor. This char- 
acteristic spectrum has remained elusive because it only 
appears at temperatures far below T c in very clean crys- 
tals, thus necessitating the development of a sensitive 
bolometric technique for making broadband microwave 
measurements. 

A simple feature of a clean d x 2_ y 2 superconductor 
is that the low energy density of states is expected 
to have the form N(e) oc e/Ao where the linear en- 
ergy dependence comes from the linear dispersion of 
the superconducting gap near its nodes in momentum- 
space and A is the energy scale of the superconduct- 
ing gap. At sufficiently low temperatures the microwave 
conductivity should be governed by the quasiparticle ex- 
citations near the nodes. The microscopic details of 
the scattering mechanism and the framework for calcu- 
lating the conductivity provide an opening where one 
might find novel physics entering into the problem. 
In the low temperature limit, the conventional frame- 
work for calculating the conductivity in the supercon- 
ducting state is a self-consistent ^-matrix approxima- 
tion (SCTMA) developed to account for impurity pair- 
breaking effects which modify both N(e) and the elec- 
tron self-energy This work showed that the ex- 
pression for the microwave conductivity for any scatter- 



ing strength takes a simple energy-averaged Drude form 
a(u),T) = ne 2 /m* ([iw + T _1 (e)] _1 ) e , where the energy 
dependence of the scattering rate r _1 (e) is determined 
by the strength of the scattering and < • • • > E denotes a 
thermal average weighted by N(e) [Q. At low T, the uni- 
tary (strong-scattering) limit result has r _1 (e) w T u /e 
while the Born limit has T _1 (e) rj TbS (to within loga- 
rithmic corrections), where the scale factors T u and Yb 
are determined by the density of impurities. The puzzling 
failure of this theory has been that neither form resembles 
the energy-independent r _1 inferred from the microwave 
measurements on YBa2Cu306.gg by Hosseini et cd. Q. 

Two major technical steps have now led to the dis- 
covery of a low temperature regime where some of the 
expectations of the simple theory seem to hold. First, 
the samples used in the experiment reported here are of 
exceptionally high crystalline perfection. Most cuprate 
materials retain a substantial level of intrinsic disorder 
due to cation cross-substitution, but the YBa2Cu306+ x 
system has a chemical stability that guarantees an ex- 
tremely low level of cation disorder (<10 -4 ). Large im- 
provements in the purity of YBa2Cu30g+x crystals were 
brought about by the advent of BaZr03 crucibles, which 
do not corrode during crystal growth ||. With a very 
high degree of atomic order in the Cu02 planes, the sec- 
ondary disorder effect of the off-plane oxygen atoms be- 
comes important. Since intercalation of oxygen into the 
CuO chains in this material (which run along the crystal 
&-axis) act as the reservoir for hole doping, probing the 
phase diagram typically involves nonstoichiometry and 
off-plane disorder. However, a few highly ordered phases 
are available in this system. Here we concentrate on fully- 
doped YBa2Cu307, with every CuO chain filled, and 
ortho-II ordered YBa2Cu30g.5o, with every other CuO 
chain filled. The perfection of the fully-doped samples 
is limited by not being able to anneal to completely full 
chains, 06. 99 being our limit. The crystalline perfection 
of the YBa2Cu306.5o crystals comes from the high pu- 
rity which in turn allows the development of very long 
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correlation lengths of the ortho-II order in the 3 crystal- 
lographic directions (£ Q = 148 A, £,=430 A, and £ c =58 A) 
P|. The measurements shown here were performed on 
mechanically detwinned single crystals of 99.995% purity 
with dimensions (a x b x c) ~ (1 x 1 x 0.010) mm 3 . 

The second key step to the results presented here is a 
measurement technique based on a bolometric method of 
detection, which provides a natural way of covering the 
microwave spectrum in more detail than is possible with 
a set of fixed-frequency experiments. The technique has 
recently proven useful for measuring resonant absorption 
in the cuprates 0. For the present work we have im- 
proved the sensitivity of the technique to the pW level 
necessary for resolving the intrinsic power absorption of 
a small single crystal, while at the same time employing 
an in-situ Ag:Au alloy reference sample of known sur- 
face resistance that calibrates the absorption absolutely. 
This is particularly important because it allows compen- 
sation for the strong frequency dependence of the mi- 
crowave field amplitude caused by standing waves in the 
microwave circuit. 

In our system, the microwave surface resistance 
R s {lu,T) is inferred from the synchronous measurement 
of the sample temperature rise as the amplitude of the 
microwave magnetic field H r f is modulated at low fre- 
quency (1 Hz). The high-T c and reference samples are 
placed in symmetric positions in a rectangular shaped 
coaxial transmission line with a broadened center con- 
ductor designed to ensure spatial magnetic field homo- 
geneity, as depicted in the inset of Fig. [l| The coaxial line 
is terminated with a short, providing a region where the 
samples experience predominantly magnetic fields. Both 
samples are mounted on thin sapphire plates with a small 
amount of silicone grease. Chip heaters mounted on the 
thermal stages allow precise calibration of the thermal 
sensitivity, while the use of a second Ag:Au sample in 
place of the high-T c sample confirms the absolute value 
and frequency independence of the calibration. All mea- 
surements use the low demagnetization factor geometry 
for a thin platelet, with H r f || b. Thus the screening cur- 
rents run along the crystal's a-axis, completing the loop 
with a short section of c-axis current which we do not at- 
tempt to correct for but know to be a small contribution. 

Figure [l] shows the frequency dependence of R s for a- 
axis currents in a crystal of YBa2Cu30g.5o- The measure- 
ments span 0.6-21 GHz, limited at low frequency where 
the small dissipation of the sample approaches the reso- 
lution limit of the experiment. The 1 GHz values for the 
rms uncertainty in R s , SR S , are about 0.2, 0.4, 0.6, and 
1.3 for T = 1.3, 2.7, 4.3, and 6.7 K respectfully. For 
the worst case of 6.7 K this translates to 8a\/a\= 0.14 
with uncertainties decreasing correspondingly at higher 
frequencies and lower temperatures. For clarity, error 
bars are omitted from the figures. At high frequencies, 
the apparatus is limited by deviations from the TEM field 
configuration in the sample region. An overall 5% uncer- 
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FIG. 1: Low temperature surface resistance data for the 
a-axis direction of ortho-II ordered YBa2CuaO6.50 obtained 
from synchronous power absorption measurements. The mi- 
crowave transmission line used for the measurements is shown 
in the inset, with the platelet-shaped high T c (black) and 
Ag:Au alloy reference (grey) samples positioned in the regions 
of homogeneous H r f- 



tainty in R s is set mainly by the Ag:Au reference-alloy 
DC resistivity measurement. 

In the limit of local electrodynamics, the relevant 
limit in our sample geometry S, the surface impedance 
Z s = R s + iX s is connected to the complex conductiv- 
ity a = o\ — i(T2 via Z s = \J\[iQLoj a. Measurements of 
both R s and X s can be used to determine a\ and 02, 
but broadband measurements of only R s present a prob- 
lem. A common solution, used in infrared measurements, 
is to measure a single quantity such as reflectance over 
a wide enough range that Kramers-Kronig and Fresnel 
equations can be used to calculate o\ and 02. Another 
solution is to fit spectra to a model of (j\ and er 2 that au- 
tomatically satisfies the Kramers-Kronig relations. How- 
ever, low frequency microwave measurements allow an- 
other option for extracting o\ (w, T) from R s (oj, T) . Deep 
in the superconducting state, both 02 and X s are largely 
controlled by the screening associated with the response 
of the superfiuid. In the low frequency limit X s = hqljX 
and (72 = (/Lt wA 2 ) _1 where A is the London penetration 
depth. So, at low temperatures a broadband measure- 
ment of R s (uj, T) together with a measurement of A(T) at 
a single frequency is sufficient to extract the conductivity 
spectrum ai(tu, T). Using cavity perturbation techniques 
with a superconducting loop-gap resonator operating at 
1.1 GHz § we measure AX S (T) = (j, uA\(T), where 
AA(T) = A(T) - A(T = 1.2 K) is the temperature de- 
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FIG. 2: The low T evolution of the quasiparticle conductiv- 
ity spectra of ortho-II YBa2Cu3Oe.50 (same symbols as Fig. 
[l]). The solid curve is a fit to the 6.7 K data with the Born- 
scattering model, but the dashed curves show that the model 
fails to capture the observed T dependence. The inset shows 
that the data obey an unusual frequency-temperature scal- 
ing o-i (w,T) = (7x(u>/[T + T ]) with T = 2.0 K. The Drude 
fit in the inset illustrates the inadequacy of the Lorentzian 
lineshape. 



pendent increase of the London penetration depth from 
its low temperature limit. Since these measurements do 
not determine A absolutely we take the T — > value to 
be Ao = 1600 A for overdoped YBa2Cu306.gg as mea- 
sured by infrared reflectance jlO) and Ao = 2600 A for 
underdoped YBa2Cu30g.5o, consistent with /xSR mea- 
surements A more accurate extraction of ai(u>,T) 
involves a correction for screening by the normal fluid 
conductivity. To do this we use a self-consistent method 
that includes contributions to 02 (w, T) both from the su- 
pcrfluid and from the screening component due to quasi- 
particles Q. At the low temperatures of this experiment, 
where there are few quasiparticles to contribute to the 
screening, this is a small correction having a maximum 
effect of 8% on <j\ for our highest temperatures and fre- 
quencies. 

Figure |2| depicts the real part of the quasiparticle con- 
ductivity extracted from the broadband R s measure- 
ments. The strong frequency dependence over intervals 
as small as 1 GHz (equivalent to a temperature of 0.05 K) 
shows that the timescale associated with the scatter- 
ing of low-energy quasiparticles in these extremely clean 
samples falls within the bandwidth of our measurement. 
The cusp-like shape of <ji(u>,T), the approximately T- 
independent low frequency limit and a tail that falls more 
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FIG. 3: Phenomenological fits with Eq. [j] (same symbols as 
Fig. 0). The inset compares the normal fluid density n n e 2 /m* 
obtained by integrating Eqn. [I] with the loss of superfluid den- 
sity inferred from 1.1 GHz AA(T) data. The dashed line has 
the slope of the open symbols, and the agreement indicates 
that the normal-fluid and superfluid spectral weight obey the 
conductivity sum rule, although a residual normal fluid term 
is implied. 

slowly than 1/ui 2 are all features expected for Born scat- 
tering. The solid curve in Fig. |^ shows a convincing fit 
to the 6.7 K data using the energy-averaged Drude form 
with r _1 (e) = TbS, indicating that the overall shape 
of the spectrum is well-described by the weak scatter- 
ing calculation with fit parameters HTb = 0.032 and 
ne 2 h/(m*A ) = 5.9 x 10 5 fi -1 m _1 . The other curves in 
Fig. y are the Born-limit predictions for the lower tem- 
peratures, using the parameters that fit the 6.7 K data. 
It is clear that this model progressively underestimates 
the spectral weight as T is reduced, and because of this 
a global fit for all temperatures produces results that 
are less satisfactory. The inset of Fig. || shows that the 
o-i (w,T) data scale as cj/(T + T ) with T =2.0 K > rather 
than scaling as uj/T as expected in the Born limit. Hence, 
the weak scattering limit captures the frequency depen- 
dence of the data, but the requirement of the SCTMA 
model that the spectral weight vanish as T — ► leads to 
disagreement. 

In order to fit to the spectra, we adopt a phe- 
nomenological form for a\ {uj, T) that captures the Born- 
lineshape features, namely the cusp-like shape and high 
frequency tail: 

ffi(w,T) = <ro/[l + (w/r)»]. (1) 

We fit directly to R s (ui,T), thereby automatically ac- 
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FIG. 4: The conductivity spectrum of a fully-doped sample 
of YBa2Cu3C>6.99 in the a-direction. The Drude fits to the 
spectra highlight the evolution from a cusp-like shape to a 
more Lorentzian lineshape with increasing temperature. 



counting for the quasiparticle screening. Fig. y shows 
the fits to individual spectra using this model where the 
parameters a and y remain relatively constant, tak- 
ing average values of 1.67(±0.05) x 10 7 Vl^ 1 m~ 1 and 
1.45 (±0.06) respectively. The parameter T varies approx- 
imately linearly in T with fit values 12.1, 19.2, 26.3, and 
35.0 x 10 9 sec _1 . Although these fits also suggest the 
unusual frequency-temperature scaling discussed previ- 
ously, enforcing it in the model reduces the agreement 
in the spectral weight comparison. The fits can be inte- 
grated in order to obtain absolutely the temperature de- 
pendent spectral weight associated with the normal fluid. 
In the inset of Fig. ^ this normal fluid spectral weight 
is compared to the spectral weight lost from the super- 
fluid as determined independently by the measurements 
of AA(T). The slopes agree to within 2%, a good indica- 
tion that the model is capturing the oscillator strength at 
higher frequencies. It is important to note that at low T 
this comparison is independent of the choice of Ao to first 
order since a x w 2R s /(^cj 2 \f l ) and A(l/A 2 ) « — AA/Aq. 
The offset apparent in the normal fluid density corre- 
sponds to a T = residual normal fluid. 

The discovery of a cusp-shaped spectrum for the low 
temperature conductivity of YBa2Cu30g.5o seems at 
odds with the Drude fits (y = 2 in Eq. Q) that were 
found to reasonably describe the spectra inferred from 
fixed-frequency microwave measurements on fully-doped 
YBa 2 Cu306.99 Q. To resolve this conflict, broadband 
conductivity spectra were obtained for the a-axis of a 
crystal of the fully doped material, shown in Fig. [|. At 



the lowest temperatures there is indeed a cusp-like spec- 
trum, similar to the one seen in the underdoped ortho-II 
sample. The reason this was not seen in earlier mea- 
surements is that the spectrum gives way to a more 
Lorentzian lineshape above 4 K, as indicated by the pro- 
gressively better fit to a Drude model with increasing 
temperature; the 5 spectra do not scale in the manner 
seen for YBa2Cu3C>6.5o- Hosseini et al. concluded that 
the scattering rate was nearly temperature independent 
below 20 K at a value of t" 1 = 5.6(±0.6) x 10 
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and here we find r 



4.4(±0.3) x 10 



10 



sec 



sec , 
1 with 



the decrease likely due to continued improvements in the 
sample purity. The integration of the Drude model cap- 
tures the oscillator strength rather well, except at the 
lowest temperatures where the fits are too poor for us 
to comment upon an extrapolated residual value. It is 
not yet clear why the non-Lorentzian lineshapes are re- 
stricted to such low temperatures in the fully doped sam- 
ple. 

To conclude, we have provided the first highly detailed 
measurement of the quasiparticle conductivity spectrum 
in the disorder-dominated regime of an extremely clean 
d-wave superconductor. This regime has been accessed 
by using well-ordered crystals with very high purity, mea- 
sured with a broadband microwave technique whose fre- 
quency range matches the very small quasiparticle scat- 
tering rate in the samples. A number of puzzles remain, 
particularly the residual normal fluid inferred from ex- 
trapolations to T = 0, a phenomenon that is seen to 
a much greater degree in other cuprates |12, The 
cause of the evolution from a cusp-like spectrum to a 
more Lorentzian lineshape in YBa2Cu306.gg is also un- 
clear and may require a more sophisticated scattering 
model that employs intermediate scattering strengths, 
neither Born nor unitary. Nevertheless, several features 
of the quasiparticle scattering dynamics are characteris- 
tic of Born-limit scattering, namely a cusp-shaped con- 
ductivity spectrum seen at low temperatures. Additional 
physics such as order parameter suppression at impurity 
sites jlill might be needed to resolve the remaining puz- 
zles, but the data presented here provide a simple start- 
ing point that is quite close to the expectation for nodal 
quasiparticles scattered weakly by impurities. 

We gratefully acknowledge useful discussions with A.J. 
Berlinsky, C. Kallin, and P.J. Hirschfeld, as well as finan- 
cial support from the Natural Science and Engineering 
Research Council of Canada and the Canadian Institute 
for Advanced Research. 



* Permanent address:Dept. of Physics, Simon Fraser Uni- 
versity, Burnaby, B.C., Canada, V5A 1S6. 

t Permanent address:Dept. of Physics and Astronomy, Mc- 
Master University, Hamilton, ON, Canada, L8S 4M1. 
[1] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 



■5 



(2000). 

[2] P. J. Hirschfeld, W. O. Putikka, and D. J. Scalapino, 
Phys. Rev. Lett. 71, 3705 (1993); Phys. Rev. B 50, 10250 
(1994). 

[3] C. T. Rieck, D. Straub, and K. Scharnberg, J. Low 
Temp. Phys. 117, 1295 (1999). 

[4] A. Hosseini et al., Phys. Rev. B 60, 1349 (1999). 

[5] A. Erb, E. Walker, and R. Flukiger, Physica C 258, 9 
(1996); R. X. Liang, D. A. Bonn, and W. N. Hardy, Phys- 
ica C 304, 105 (1998). 

[6] R. X. Liang, D. A. Bonn, and W. N. Hardy, Physica C 
336, 57 (2000). 

[7] Y. Matsuda et al., Phys. Rev. B 49, 4380 (1994); O. K. C. 



Tsui et al., Phys. Rev. Lett. 73, 724 (1994); M. Gaifullin 
et al., Phys. Rev. Lett. 83, 3928 (1999). 
[8] I. Kosztin and A. J. Leggett, Phys. Rev. Lett. 79, 135 
(1997). 

[9] W. N. Hardy et al., Phys. Rev. Lett. 70, 3999 (1993). 
[10] D. N. Basov et al., Phys. Rev. Lett. 74, 598 (1995). 
[11] R. I. Miller et al. (Unpublished). 
[12] D. M. Broun et al. (Unpublished). 
[13] J. Corson et al., Phys. Rev. Lett. 85, 2569 (2000). 
[14] M. H. Hettler and P. J. Hirschfeld, Phys. Rev. B 61, 
11313 (2000). 



